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significant  importance  in  vibration  COLLISION  PROCESSES  aNu 

:esses  at  thermal  velocities* 


SECTION  I 


INTRODUCTION 

The  purpose  of  this  contract  is  to  develop  a theory  that 
describes  the  vibrational-vibrational  collision  processes  within  the 
bending  mode  of  CO2  and  the  cross  collisions  between  the  bending  mode 
and  the  symmetric  stretch  mode.  These  processes  play  a vital  role  in 
the  deactivation  of  the  lower  state  of  the  CO2  laser  transition.  The 
technical  problems  associated  with  the  problem  are  simply  stated.  First, 
CO2  is  a linear  triatomic  molecule,  and  it  is  unclear  whether  e.xistent 
theories,  which  are  restricted  to  diatomic  systems,  will  work.  Second, 
since  CO2  is  linear,  the  bending  mode  is  degenerate,  which  leads  to  an 
ambiguity  between  vibrational  and  rotational  motion  (called  "vibrational 
angular  momentum"),  which  gives  rise  to  a doubling  of  the  bending  mode. 
Because  of  the  uncertainty  with  how  to  proceed,  a Bom  thec^y  is  used, 
which  is  relatively  simple  conceptually  and  which  can,  in  principle, 
lead  to  the  information  needed  to  rationally  attack  the  problem. 
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TECHNIQUE 

In  this  section,  a simple  technique  is  presented  for  calculating 
the  matrix  elements  of  the  intermolecular  potential  between  colliding 
molecules  undergoing  vibrational  energy  transfer  or  conversion  pro- 
cesses. The  theory  of  vibrational  relaxation  processes  in  diatomic 
or  polyatomic  molecular  gases  (Refs.  1 and  2)  depends  critically  on 
these  matrix  elements.  To  date,  most  calculations  of  (molecular) 
inelastic  collision  cross  sections  have  employed  a set  of  questionable 
approximations  in  computing  the  relevant  matrix  elements  (Refs.  3 
through  S) . While  these  approximations  greatly  simplify  the  calcula- 
tions, they  lack  sound  justification  and  in  some  cases  are  not  valid. 

For  collisions  between  two  molecular  systems,  most  calculations  consider 
only  (1)  collinear  collisions  and  (2)  interactions  between  nearest 
neighbors.  In  addition,  the  colliding  potential  is  approximated  by 
a simple  repulsive  exponential . 

The  first  two  approximations  are  kinematic,  and  their  effect  is 
to  eliminate  all  of  the  angular  factors  that  would  otherwise  arise  in 
the  collision  process.  To  simulate  the  full  set  of  possible  initial 
collision  orientations,  a steric  factor  Pq  = 1/5  for  each  molecule  is 
introduced  into  the  transition  amplitude.  The  second  approximation 
is  justified  on  the  grounds  that  it  is  the  short-range  portion  of  the 
collision  potential  that  dominates  the  dynamics  of  vibrational  transfer 
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and  conversion  processes.  Unfortunately,  the  steric  factor  tends  to 
vary  with  the  colliding  molecular  species,  and  there  does  not  exist  a 
reliable  quantitative  theory  for  calculating  this  factor  for  collisions 
involving  nonlinear  polyatomic  molecules  (Ref.  6) . The  second  approxi- 
mation is  clearly  invalid  for  initial  collision  orientations  in  which 
one  atom  of  one  molecule  is  very  close  to  a number  of  atoms  in  the 
second  molecule.  Furthermore,  one  cannot  investigate  the  possibility 
of  rotational  transfer  or  conversion  occiirring  simultaneously  with 
vibrational  processes.  Finally,  with  regard  to  these  kinematic  approxi- 
mations, it  should  be  noted  that  Mies  (Ref.  7)  has  carried  out  a 
first-principle's  calculation  using  the  full-angle-dependent  potential 
for  the  excitation  of  hydrogen  molecules  by  helium  atoms.  His  results 
have  raised  grave  doubts  about  the  validity  of  items  (1)  and  (2)  above. 
Assume  that  the  potential  can  be  approximated  by 

U = lU^  At  - (1) 

Here,  U.  .is  the  interaction  between  atom  i in  polyatomic  molecule  A 

and  atom  j in  polyatomic  molecule  B,  r is  the  vector  between  the  molecular 

center  of  masses,  and  is  the  position  vector  of  atom  i(j) 

t J 

from  the  center  of  mass  of  molecule  A(B).  Equation  (1)  assumes  that 


the  electronic  degrees  of  freedom  are  not  excited  by  the  collision  and 

in  addition,  that  s » r - is  the  only  vector  direction 

V j 

that  is  significant.  Thus,  one  neglects  (orientational)  dipole-dipole 
forces,  etc.,  which  are  believed  to  dominate  only  the  long-range  portion 
of  the  intermolecular  potential  (Ref.  8) . This  part  of  the  intermolecular 
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force  gives  rise  only  to  rotational  processes  and  plays  no  role  in 
vibrational  relaxation. 


For  situations  in  which  the  Bom  approximation  applies,  the  matrix 
elements  of  equation  (2)  are  required  between  the  states  |^;a,3>  and 

Here  is  the  relative  center-of-mass  wave  vector  before 

(after)  the  collision,  a(a')  is  the  initial  (final)  rotation-vibration 
state  of  molecule  A,  and  3(3')  is  the  same,  but  for  molecule  B.  The 

relevant  matrix  element  is  then  j 

<l^';a',3'lyl^;a,3>  = .1  (2) 

i,d 

where -ftq  is  the  momentum  transfer  vector  and  = <a' | exp(-iq*R^^^^  la>  / 

is  the  form  factor  of  atom  i between  molecular  states  |a>  and  |a'>  at 

wave  vector  q.  These  quantities,  which  also  appear  in  neutron  scattering, 

depend  only  on  the  molecular  structure,  the  initial  and  final  rotation-  ' 

vibration  states,  and  the  wave  vector  q.  They  are  independent  of  the 

colliding  potential  and  for  some  transitions  may  be  obtainable  from 

neutron  scattering  experiments  (Ref,  9).  The  factors  .j 

Fourier  transform  of  the  atom-atom  interaction.  This  approach  separates  J 

the  internal  dynamics  of  the  molecule  (i.e.,  the  form  factors  C-<l)) 

from  the  center-of-mass  motion  (i.e.,  the  U.  .(q))  and  is  formally 

J 1 

similar  to  the  approach  taken  for  inelastic  collisions  between  fast 

electrons  and  atoms  (Ref.  10).  Finally,  it  should  be  noted  that  all  of  j 

the  approximations  mentioned  above  have  been  avoided. 

I 

For  most  molecular  systems  one  cannot  use  the  Born  approximation  ] 

i 

for  collisions  occurring  at  thermal  velocities  and  one  must  resort  to 
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other  techniques.  One  approach  is  the  Born  approximation  to  the  reactance 
matrix.  This  technique  was  used  by  Burke  and  Seaton  (Ref.  11)  to  study 
e-H  collision  cross  sections.  Their  efforts  were  successful  and  can  be 
adapted  readily  to  our  problem. 

In  general,  one  expects  that  many  partial  waves  will  be  required 
to  adequately  describe  the  cioss  section  for  inelastic  processes.  Con- 
sequently, one  should  avoid  a partial  wave  analysis  and  use  the  expression 

1 - ■i(2mk/fi^)<k' ,a' ,3' |y|  k,o,6> 


for  the  transmission  matrix. 
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SECTION  III 


DISCUSSION 


As  discussed  in  Section  II,  the  collision  matrix  elements  decompose 
into  individual  atom-atom  contributions  times  a term  describing  the 
center-of-mass  motion.  This  leads  to  a model  of  collisions  that  is 
called  the  "form  factors"  model  following  the  usage  of  this  model  in 
nuclear  theory.  The  basic  idea  is  to  regard  each  pair  of  components 
(e.g.,  0 and  C)  to  be  the  same  independent  of  the  molecule  in  which  it 
is  situated.  In  practice,  the  atmosphere  has  only  four  interesting 
atoms  (H,  0,  C,  and  N)  and  the  CO2  laser  only  one  more  (He);  thus,  one 
can  describe  all  collisions  of  interest  with  a relatively  limited 
number  of  parameters.  In  the  case  of  the  collision  of  CO2  molecules, 
there  are  only  three  collision  potentials  (0-0,  0-C,  C-C)  that  can  be 
determined  by  comparison  with  the  collisions  of  CO  with  itself--a 
problem  that  is  fairly  well  developed  both  theoretically  and  e.xperimentally. 

In  the  specific  approach,  a relatively  simple  problem  was  chosen, 
the  collision  of  He  with  N2,  because  of  its  significance  to  the  CO2 
laser  problem  and  because  it  involves  one  of  the  simplest  of  the 
homonuclear  diatomics.  Furthermore,  it  is  somewhat  analogous  to  colli- 
sions involving  hydrogen,  which  is  a system  about  which  a great  deal  is 
known.  This  work  was  initially  developed  into  a short  paper  whose 
material  is  included  in  this  report.  Several  purposes  were  served 
in  the  study.  It  showed  that  the  results  of  the  calculation  were 
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sensible  and  that  calculations  were  possible.  More  subtly  it  demonstrated 
that  quantum  mechanical  effects  were  important. 

The  investigation  given  in  the  paper  involved  relative  cross  sections 
which  are  used  to  determine  the  range  but  not  the  magnitude  of  the 
potential.  For  that,  one  must  be  able  to  calculate  the  absolute 
cross  sections.  To  do  this,  it  is  necessary  to  unitarize  the  calculation, 
i.e.,  to  correct  for  the  fact  that  probability  conservation  is  violated 
in  the  Born  approximation.  Many  unitarization  schemes  exist  and  are 
tested  for  collisions  with  hydrogen.  The  scheme  outlined  in  Section  II 
was  tried  and  proved  to  be  successful  for  e~  and  He  collisions.  How- 
ever, an  unresolved  difficulty  was  encountered.  The  unitarization  pro- 
cedure involves  the  computation  of  all  nonunitarized  cross  sections 
before  the  unitarized  values  can  be  found.  As  is  seen  in  Section  IV, 
there  are  significant  contributions  from  large  values  of  angular  momentum 
transfer  (A*/"  = 12)  and  therefore  intricate  dependences  on  the  m values 
(i.e.,  angular  dependences).  These  numerical  problems  become  worse  as 
J increases,  and  there  are  enormous  problems  in  finding  an  appropriate 
cutoff  point.  At  the  end  of  the  contract  year,  this  problem  had  not 
been  solved.  However,  unitarization  may  not  be  as  important  as  first 
believed,  and  the  difficulties  may  not  be  fatal  to  the  scheme. 
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SECTION  IV 


CALCULATIONS 

In  this  section,  the  method  described  above  is  applied  to  calculate 
the  scattering  cross  sections  for  vibrational,  rotational,  translational 
collision  rates  between  He  and  N2.  In  particular,  excitation  from  the 
vibrational  ground  state  (v  = 0)  to  the  v = 1 state  is  considered  and  the 
following  three  features  are  examined:  (1)  The  effect  of  the  initial 
rotational  state,  (2)  the  importance  of  angular  momentum  transfer  during 
a vibrational  encounter,  and  (3)  the  effects  of  molecular  orientation 
Ci.e.,  the  Mj  quantum  number)  in  these  collisions. 

For  the  purposes  of  illustration,  the  inelastic  scattering  of  He 
with  N2  is  treated  within  the  first  Bom  approximation.  As  inelastic 
vibrational  processes  generally  involve  the  short-range  portion  of  the 
collision  potential,  which  is  quite  large,  the  first  Bom  approximation 
will  not  yield  collision  cross  sections  whose  absolute  magnitudes  are 
reliable.  However,  it  does  yield  accurate  results  for  the  relative 
sites  of  such  cross  sections  for  different  inelastic  processes.  For 
example,  in  charge  transfer  collisions  between  alpha  particles  and 
hydrogen,  the  ratios  of  the  inelastic  cross  sections  to  different  final 
states  remain  unchanged  in  going  from  the  Bom  approximation  to  much 
more  sophisticated  approaches  although  the  magnitude  of  the  various 
cross  sections  varies  by  several  factors  of  10  (Refs.  12  and  15). 
Furthermore,  experimental  results  for  these  ratios  are  within  10°6 


agreement  with  those  predicted  by  the  Born  approximation.* 

For  our  present  purposes,  these  ratios  are  sufficient,  and  our  results 
are  qualitatively  reliable. 

Using  the  techniques  described  in  Section  II,  it  can  be  assumed 
for  the  purposes  of  vibrational  collisions  that  the  short-range  piece 
of  the  He-N2  interaction  can  be  written  as 

- ti)  * 1/2 (r  - ^2)  (4) 


where  r is  the  vector  distance  between  the  He  atom  and  the  center  of 
mass  of  the  N2  molecule.  The  vectors  Si  and  S2  locate  the  position  of 
the  nitrogen  atom  in  question  relative  to  the  molecular  center  of  mass. 
The  quantity  i'l (r  - Si){£/2(r  - S2)]  is  the  interaction  potential  between 
atom  1[2]  and  the  He  atom.  For  homonuclear  diatomic  molecules  l/i  ■ i/2» 
and  for  simplicity  a Yukawa  type  interaction  with  a range  of  0.1  nm 
has  been  assvimed. 

Within  the  first  Bora  approximation,  the  inelastic  scattering 
cross  section  between  the  (initial)  molecular  state  and  the 

(final)  molecular  state  is  (the  He  atom  is  assumed  to  remain 

in  its  ground  state) 


J V 


V 


dgt/i^d)  \y 


e/  t 


. (g) 


(5) 


where  p is  the  reduced  mass  of  the  He-N2  system,  and  O'i(g)  is  the  Fourier 
transform  of  the  Yukawa  potential.  For  a homonuclear  diatomic  molecule 
such  as  N2,  the  form  factor 

J J 
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Yale  data,  unpublished. 


'yw«,*vV«  ■'*)  • 2<vV«_,'|cos(ig-i5)|vJ«_,> 
j j 

where  R is  the  instantaneous  vector  distance  between  the  two  nitrogen 

atoms  (in  the  N2  molecule) . If  the  axis  of  quantization  is  taken 

parallel  to  the  momentum  transfer  vector  (g) , then  the  selection  rule 

M,  = Wt'  follows.  A similar  situation  arises  in  neutron  scattering 
«/  «/ 

(Ref.  9).  Finally,  note  that  symmetry  demands  that  both  J and  J'  be 
even. 
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SECTION  V 

RESULTS 


For  convenience,  this  discussion  is  confined  to  vibrational  excita- 
tion (de- excitation)  from  the  ground  (first-excited)  to  the  first- 
excited  (ground)  vibrational  state.  In  problems  involving  atmospheric 
propagation  or  in  typical  laser  systems  one  generally  finds  that  a 
large  number  of  different  rotational  states  are  occupied.  This  is 
especially  true  for  molecules  such  as  CO2  or  N2.  Consequently,  one 
must  investigate  the  effect  of  the  initial  angular  momentiM  state  on 
the  collision  process.  One  can  readily  convince  oneself  that  the  best 
measure  of  this  feature  is  to  examine  V-»-T  collision  process,  i.e., 
inelastic  events  in  which  i^J  = 0 (but  with  Av  = 1 , of  course) . In 
Figure  1 the  collision  cross  section  for  = 0 to  = 10  is  plotted  in 
arbitrary  units  for  an  incoming  velocity  of  v * 1.14  v^,  where  is 
the  threshold  velocity  for  such  inelastic  processes.  For  J > A,  the 
total  inelastic  cross  sections  scale  with  the  multiplicity  of  the  state, 
which  implies  that  the  probability  for  inelastic  scattering  (with 
AJ’  = 0)  is  independent  of  J for  the  higher  rotational  states.  This 
feature  is  unaffected  by  increasing  the  initial  velocity.  In  addition, 
by  V » 1.8  V , the  J - 2 state  grows  until  it  scales  with  its  multi- 
plicity,  and  for  very  high  velocities  (v  « 4v^)  the  ground  rotational 
states  also  do  the  same.  In  addition,  this  feature  is  further  rein- 
forced by  the  fact  that,  at  a distance  from  threshold, 


l.Or 


OS 


V l.l4Vj 


I 


8 


9 


Figure  1.  Plot  of  normalized  calculated  cross 
sections  (in  arbitrary  units)  as  a 
function  of  J. 
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J- 


(7) 


for  all  states  through  J = 12.  In  equation  (7),  a^^,(a^,^)  is  the 
total  inelastic  scattering  cross  section  for  processes  in  which  Av  = 1 
and  the  rotational  quantum  number  changes  by  means  of  J ■*  J) . 

(The  small  deviation  from  unity  is  due  to  small  changes  in  the  available 
phase  space  for  collision.) 

Next,  the  importance  of  angular  momentum  in  a vibrational-roational- 
translational  (V-*11->T)  collision  is  exa;ained.  Processes  in  which  ^ 0 
represent  additional  channels  through  which  the  vibrational  mode  can 
be  excited  or  de-excited  and  as  a result  may  play  an  important  role 
in  gas  kinetics.  This  will  occur  if  the  cross  sections  for  Ae/'  ^ 0 
processes  are  of  the  same  order  of  magnitude  as  the  V « 0 channel.  In 
Figure  2 the  quantities 


'^i.2-K).cr^''^ 

T? 

I i^i  -(v) 


(8) 


12 


{ 


^ 5 io  S 20  25  30  35  40  45 

Velocity  (in  10^  cm/ sec) 

I ' ’ 

! 

Figure  2.  Plot  of  normalized  calculated  cross  sections 
(in  arbitrary  units)  for  de-excitation  of  N2 
by  He  as  a function  of  velocity. 

9 

are  plotted  for  de-excitation  of  N2  by  He  as  a function  of  v,  the 
incident  velocity.  Note  that  large  quantities  of  angular  momentum 
are  readily  transferred.  In  particular,  Ae7  = 0 processes  amount  to 

i 

less  than  10%  of  the  possible  processes  for  thermal  velocities  of  1200®K  i 

or  greater.  In  addition,  for  thermal  velocities  of  SOO^K,  processes  \ 

with  = 12  are  of  greater  importance  than  those  in  which  A«7  = 0. 

Furthermore,  for  a cold  gas  (105®K) , processes  with  iW  = 14  are  as 
significant  as  those  with  A«7  = 0.  Thus,  angular  momentum  transfer  plays 
an  extremely  important  role  in  vibrational  excitation.  Note  that, 
except  for  the  c/  = 2 -►  «/  = 0 processes,  these  effects  are  not  due  to 
multiplicity  of  the  final  state  but  arise  from  the  selection  rule 
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tiMj  = 0 state.  Therefore  only  one  Mj  (i.e.,  Mj  = 0)  state  can  par- 
ticipate, which  gives  rise  to  an  anomalously  small  cross  section.  In 
Figure  3 a similar  plot  has  been  made  for  vibrational  excitation  includ- 
ing states  through  J = 12.  For  convenience,  the  relative  cross  sections 
have  been  plotted  as  a fiinction  of  v/v^,  where  is  the  threshold 
velocity  for  excitation  of  N2  by  He  to  the  first  excited  vibrational 
state.  The  situation  is  essentially  identical,  although  a much  greater 
range  of  impact  velocities  is  displayed. 


Figure  3.  Plot  of  normalized  calculated  cross  sections 
(in  arbitrary  units)  for  excitation  of  N2  by 
He  as  a function  of  velocity. 


An  examination  of  Figure  3 reveals  that  as  the  incident  velocity  increases, 
the  V = 0 channel  tends  to  dominate  and  the  probability  of  creating 
large  quantities  of  angular  momentum  tends  to  decrease.  However,  even 
at  an  incident  velocity  of  4v  , the  AJ  = 0 channel  accounts  for  only 
half  the  possible  processes,  so  V-^R-*-T  processes  are  always  of  importance. 
These  features,  i.e.,  large  Ai7,  are  quantum  mechanical  effects,  which 


would  not  appear  in  semiclassical  theories.  This  can  be  seen  from  the 
following  semiclassical  argument.  The  torque  ? exerted  on  the  N2  mole- 


cule is 

f . I . J . f (9) 

where  r is  the  impact  parameter  and  must  lie  within  the  range  a of  the  force. 

a 

The  change  in  the  rotational  motion  is  approximately  r 1 

torque,  i.e.,  j , 

SJ  = Fa^/r  (10) 

In  equation  (7),  v is  the  average  of  the  initial  and  final  velocities. 

Thus,  on  the  basis  of  semiclassical  considerations  one  expects  that  ' 

tj  should  be  small  and  that  the  cross  sections  with  J ^ J'  should 

fall  off  rapidly  with  increasing  velocity.  These  features  tend  to 

appear  only  at  high  incident  velocities  (v  'v  4v^) , but  note  that  even 

at  V = 5.5  quantum  effects  still  remain. 

Finally,  the  effect  of  molecular  orientation  on  V-*-T  processes  is 
examined.  For  large  values  of  J , the  angular  momentum  behaves  like  a ; 

classical  variable  (with  regard  to  the  orientation  of  the  molecule)  and 
one  can  make  a correspondence  between  the  Mj  state  and  a classical 
angle  9 h cos' ^ , which  specifies  the  orientation  of  the  molecule. 

1 

Now,  for  situations  in  which  the  incident  velocity  is  slightly  greater  | 

I 

than  the  threshold  velocity,  the  momentum  transfer  vector  is  nearly  i 

parallel  to  the  initial  wave  vector,  i.e.,  the  incident  velocity.  | 

j 

Since  the  momentiom  transfer  vector  defines  the  axis  of  quantization,  j 
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it  follows  that  states  with  Mj  = 0 correspond  to  head-on  or  collinear 
collisions  (see  inset  in  Fig.  4),  whereas  states  in  which  Mj  - ±J 
correspond  to  the  case  in  which  the  N2  molecule  is  normal  to  the  incident 
velocity.  In  Figure  4,  the  inelastic  collision  cross  section  is  plotted 
in  arbitrary  units  for  the  V-<-T  process  with  J - \2  and  Av  = 1,  as  a 
function  of  Mj  (or  6) . The  large  central  peak  corresponds  to  the  fact 
that  head-on  collisions  tend  to  dominate  V-*-T  rates.  Furthermore,  the 
configuration  in  which  the  molecule  is  normal  to  the  incident  velocity 
(i.e.,  Mj  - ±J)  does  not  participate  in  the  collision.  These  features 
are  expected  within  a classical  model.  However,  the  cross  section  does 
not  fall  off  smoothly  as  \Mj\  varies  from  zero  to  J , but  instead  dis- 
plays considerable  structure  with  intermediate  peaks.  Again,  these 
features  arise  from  quantum  effects. 


t 

b V-Vf 
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Figure  4.  Plot  of  normalized  calculated  cross  sections 


(in  arbitrary  units)  as  a function  of  the  M 
quantum  number. 
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In  conclusion,  changes  in  the  molecular  rotation  are  of  significant 


importance  in  vibration  collision  processes  and  dominate  deactivation 
processes  at  thermal  velocities. 
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SECTION  VI 

FIFTH  WINTER  COLLOQUIUM 

As  required  by  the  Statement  of  Work,  the  Fifth  Winter  Conference 
on  Quantum  Electronics  was  held  from  February  1-6,  1975,  at  Snowmass, 

Colorado.  The  emphasis  was  on  electronic  tremsitions , chemical  and 

collisional  pumping  of  electron  transitions,  and  short  wavelength  | 

optics.  Particular  emphasis  was  paid  to  molecular  physics.  E-collision 
pumping,  excitation  mechanics,  short  wavelength  lasers,  and  short 

wavelength  resonators.  • • , 
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